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ABSTRACT 


An analysis yielding six-degree-of -freedom 
equations of motion is presented for predicting the 
dynamic behavior of a general parachute-payload 
system. The parachute canopy and associated air- 
mass are approximated as a rigid body, and separate 
equations of motion are derived for the canopy and 
payload subject to the constraint of the risers and 
suspension lines. The analysis determines the forces 
and the response of various riser and suspension- line 
geometries subjected to large displacements, under 
the assumption that these lines are linearly elastic. 
The equations are readily adaptable to computer 
solutions and should be of interest in analyzing the 
dynamic performance of lifting-parachute payload 
systems. 


A MATHEMATICAL MODEL FOR CALCULATING THE FLIGHT DYNAMICS 


OF A GENERAL PARACHUTE -PAYLOAD SYSTEM 

By Joe D. Gamble 
Manned Spacecraft Center 


SUMMARY 


An analysis yielding six-degree-of-freedom equations of motion is presented for 
predicting the dynamic behavior of a general parachute -payload system. The parachute 
canopy and associated airmass were approximated as a rigid body, and separate equa- 
tions of motion were derived for the canopy and payload subject to the constraint of the 
risers and suspension lines. Equations were developed for solving for the force distri- 
bution in space frames composed of various riser and suspension- line geometries under 
the assumption that these lines are linearly elastic. The Newton-Raphson iterative 
solution was used to solve equations for large displacements, and the solution is readily 
adaptable to computer solutions. The model should be of interest for analyses of the 
performance of lifting-type parachutes, particularly when the six-degree-of-freedom 
mass and aerodynamic properties of these parachutes have been better defined. 


INTRODUCTION 


The proposed use of lifting parachutes for future land landing of manned space- 
craft will impose a greater requirement on computer simulations of the flight dynamics 
of the parachute -pay load system. Dynamic analysis requirements for the present 
water- landing systems are primarily to determine that the system is stable and that the 
oscillation limits during descent are acceptable for water-entry of the spacecraft. 

Although full-scale simulations have primarily been used to evaluate the dynamic 
performance of manned parachute systems in the past, a few attempts have been made 
to theoretically predict the dynamics. According to Ludwig and Heins (ref. 1), the first 
attempt to consider dynamic stability of a parachute-payload system was performed by 
Brodetsky (ref. 2), but few analyses followed in the interim period to 1960. Since then, 
considerable work has been done at the University of Minnesota under the direction of 
H. G. Heinrich (refs. 3 and 4), by Ludwig and Heins (refs. 1 and 5), and by personnel 
of a NASA contractor (refs. 6 and 7). 

Most of the dynamic analyses have considered the parachute-payload system to be 
a single rigid body and have linearized the equations of motion, which limits the oscilla- 
tion amplitudes to small values. Ludwig and Heins employed computer solutions to 


solve the six-degree-of -freedom motion, using the nonlinear equations of motion; but 
their analysis also assumed a single rigid-body system. 

Proposals have been made for using a lifting parachute with aerodynamic refer- 
ence lengths on the order of 100 feet which would require suspension lines of 100 feet 
or more in length. Since the suspension lines used for most parachute -payload systems 
have a relatively low spring constant, it is possible for a suspension line to lengthen 
several inches without appreciably increasing the forces or moments on the parachute. 
The suspension-line geometry is usually such that the vertical component of the 
suspension- line length is large compared to the lateral component. This geometry 
permits large lateral displacements between the payload and parachute with only small 
changes in the suspension- line lengths and thus only small changes in the forces and 
moments acting on the parachute. Since lifting parachutes are designed for lateral ma- 
neuvers, large lateral displacements between the parachute and payload are likely. 
These displacements result because the high inertia of the parachute (which must in- 
clude the apparent inertia effects of the surrounding and included airmass) will require 
large rolling moments to produce a roll attitude change during yaw maneuvers, and the 
momentum of the payload will be such as to produce a relative lateral displacement 
during the yaw maneuver. The relative displacements will have an effect on the accu- 
racy of studies which assume the parachute and payload are a single rigid body, since 
the inertia matrix will be significantly affected by such displacements. This suggests 
the need for considering the parachute and payload as separate, elastically coupled sys- 
tems when analyzing the dynamics of a lifting-parachute system. 

Another advantage of such an approach is that the aerodynamic damping of the 
parachute can be considered separately. This can result in a more accurate analysis, 
since the instantaneous aerodynamic angles of the parachute are calculated; thus, the 
effect of the damping terms C and C, are much less than when the oscillation 

m q P 

center is the combined center of mass of the parachute -payload system. With the single 
rigid-body approach, it is also impossible to account for the effect on aerodynamic 
damping of coupling of the angular rates caused by relative rotation between the two 
bodies. The two -body approach also allows the input of wind gusts acting at the actual 
parachute altitude instead of at the combined center of mass. This could be important 
in landing studies. 

Neustadt, et al. , (refs. 6 and 7) have extended the parachute analysis to include 
relative motion of the parachute and payload, with the elastic effect of the risers and 
suspension lines included. However, the analysis is restricted to planar motion and 
only considers very simplified riser geometry. 

The need obviously exists for an analysis which will predict the dynamics of a 
general parachute -payload system. This paper develops a mathematical model for an- 
alyzing the dynamics of such a system. The technique used assumes the parachute and 
payload to be two independent rigid bodies subject to the constraints of the suspension 
lines. 


The forces in the individual suspension lines and risers are calculated by solving 
the indeterminant space frames composed of these lines. These forces are then used to 
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calculate the forces and moments acting on the parachute and payload. The six-degree- 
of-freedom equations of motion of the parachute and payload are then solved using these 
forces and the aerodynamic forces. 

The number of individual space frames included in an analysis is limited only by 
the storage capacity and time requirements for a particular computer. A FORTRAN IV 
computer program developed at the Manned Spacecraft Center (MSC) provides for four 
space frames of up to 10 members each. 

The author wishes to acknowledge the assistance of William E. Thomas, MSC 
Structures and Mechanics Division, and Theodore F. Hughes, ITT/Federal Electric 
Corporation, in the preparation of this paper. 


SYMBOLS 


The following symbols are referenced to the lifting-parachute (parawing) system 
as shown in figures 1 to 5. In general, aerodynamic coefficients are functions of both 
angle of attack and angle of sideslip and are input to the computer in tabular form. If 
control inputs are available, they can be accounted for in the tables. 


[A] direction cosine matrix 

T 

[A] transpose of direction cosine matrix 

C £ aerodynamic rolling moment coefficient 


P 



oo 


dC„ 


qd 

2V„ 


9C„ 


rd 

2V 



aerodynamic pitching moment coefficient 
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m 


ac 


m 

Pd 

2V„ 


m 


0C 


m 

qd 

2 V 


n 


aerodynamic yawing moment coefficient 


n 


ac 


n 


Pd 
2 V 


ac 


n 


rd 

2V„ 


aerodynamic axial force coefficient 
aerodynamic side force coefficient 


ac 


Pd 

2V_ 


ac 


rd 

W~ 


aerodynamic normal force coefficient 


ac 


Pd 

2V„ 


ac 

z_ 

_qd_ 

2V_ 


4 


d 

M 

dx.^dz. 
dx o , dy Q , dz Q 
[F] 

F F F 
1 ’ 2 ’ 3 

Fx., Fy., Fz. 
l’ J i’ x 

Fx , Fy , Fz 
o’ J o’ o 


F A , F A , F A 
x y ’ z 

F C , F C , F C 
x y z 

G , G , G 
x’ y’ z 

I ,1 ,1 

xx yy zz 


I ,1 ,1 

xy xz’ yz 


V ■'2’ k 2 


i i £ 

l’ j l’ 1 


i i £ 
2’ J 2’ 2 


[K] 


aerodynamic reference length 

column matrix composed of elements dx^, dy., dz^ 
displacement of point i with respect to initial unstressed position 
displacement of point o with respect to initial unstressed position 
force matrix 

functional equations for Fx q , Fy Q , Fz q 

force components at point i required to hold frame in equilibrium 
virtual forces applied at frame junction 

body aerodynamic forces 

body cable forces 

gravity components along body axes 
moments of inertia about body axis system 
products of inertia about body axis system 
arbitrary frame member 

fixed points on parachute with respect to parachute axis system 

unit vector in payload axis system 

unit vector in parachute axis system 
matrix of the spring constants 
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K. 
x 

k 

L. 
x 

N 

M 



M C , M C ,M C 
x ’ y ’ z 


m 

n 

o 


spring constant of frame member i 

current step in the iteration process 
current length of frame member i 

unstressed length of frame member i 

row matrix composed of unstressed lengths of frame members 
mass 

body aerodynamic moments 

body cable moments 
total members in frame 

number of members from payload to junction point o 

space frame confluence point 

column matrix composed of the element P. 


p. axial force in member i 

x 

p, q , r body axis components of roll, pitch, and yaw inertial angular 

velocities = u> ,oj ,<x> 

X I y I Z I 



r ix, 2’ r iy, 2’ r iz, 2 


r. 

x 


free stream dynamic pressure 

inertial position vector of payload center of gravity 

inertial position vector of parachute center gravity 

scalar components of a vector from parachute to attach points on 
payload in parachute axis system 

vector from parachute to payload attach point 
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aerodynamic reference area 
time 


i m ■ ■■■■ mi hi in inn ii i ■mm 


i i hi ill ■■in mi 


i 


v x,i’ v y,i’ v z,i 

x, y, z velocity components of payload center of gravity along 
payload axes 

V x,2’ V y,2’ V z,2 

x, y, z velocity components of parachute center of gravity along 
parachute axes 

Voo 

free stream velocity of center of gravity of payload or parachute 

? 1 

inertial velocity vector of payload center of gravity 

V 2 

inertial velocity vector of parachute center of gravity 

X., Y., Z. 

l’ l l 

coordinates of attach point i 

X. , Y. ,Z. 

l o *0 *0 

original X, Y, Z coordinates of attach point i 

X , Y , Z 
°0 °0 °0 

original X, Y, Z coordinates of point o 

X Bj* Z Bj 

inertial velocity of body center of gravity along body axes 

X Bj* ^Bj’ fe Bj 

inertial acceleration of body center of gravity along body axes 

a t'*t’ y i 

X, Y, Z direction cosines of frame member i in parachute axis 
system 

[r] 

3 by 3 transformation matrix for transforming from axis system 
of body 1 to axis system of body 2 

[r] T 

transpose of [r] 

r. . 
i»] 

elements of matrix r 

AX. 

i 

X. - X - dx 
i o Q o 

AY. 

l 

Y. - Y - dy 
i o Q *o 

AZ. 

l 

Z. - Z - dz 
i o n o 


A0 


total angle change of frame member i from unstressed position 


7 



6 


correction to previous estimate in Newton-Raphson iterative 
solution 


6 correction to previous estimate of frame member i 

i 




angles between frame member i and parachute X-, Y-, and 
Z-axes 


£ } 0. original angles between frame member i in unstressed position 

*0 *0 *0 and parachute X-, Y-, and Z-axes 


p.) r., c. 
i i 


/ COS £. 

\ 1 


cos 


V /cos l //. - cos 0. V /cos 0. - COS 0. 

1 o/ V 1 V V < 


p ix, l ,P iy, l’ p iz, 1 


X, Y, Z coordinates of payload attach point of riser i in payload 
axes system 


p ix,2’ p iy,2’ p iz,2 


X, Y, Z coordinates of parachute attach point of suspension line i 
in parachute axes system 


(2 2 
\p + r + ct J 

p i 

0 

0 - • 
i] 

0 , 0,0 


0 , 0,0 


w x, 1’ w y, l’ W z,l 


w x,2 >a, y, 2 >W z,2 


w. 


column matrix composed of the elements p, r, and a (eq. (21)) 

vector from center of gravity to riser or suspension-line attach 
point 

3 by 3 matrix composed of elements of the Jacobian 
elements of the Jacobian 

relative Euler angles of body 2 with respect to body 1, applied in 
the rotation 0,0,0 

time rate of change of 0,0,0 

angular velocity components of payload about payload axes 
angular velocity components of parachute about parachute axes 
angular velocity vector of payload 


oo 


2 


angular velocity vector of parachute 



1’ 



1* 


oo 

Z, 


1 


angular velocity component of payload transformed to parachute 
axes 
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inertial angular velocity components about body axes = p, q, r 




time rate of change of velocity components about body axes 


ANALYSIS 


Equations of Motion 

The basic approach to this analysis is to assume that the parachute can be ap- 
proximated as a rigid body, and then to calculate the dynamics of the parachute and 
payload separately, subject to the constraint of connecting lines. Six-degree-of- 
freedom equations of motion can be written for each body if the forces and moments 
caused by aerodynamics as well as the forces and moments in the connecting lines are 
known. Integration of these equations will yield the motion of each body. 

The three translational equations of motion for the center of mass of a rigid body 

are 


f \ 
x 




B t 


( 


F A + F C 
x x 


F v A + F v C M 


f \ 


Zb i 

V V 


F A + F C 
z z 


/ 


G > 

y < 


z 

\ J 


0 -r q 
r 0 -p 
-q p 0 


/ \ 
x b t 


V > 


J]3 I 

V L J 


(i) 


The aerodynamic forces are determined by 
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The rotational equation of motion about the body axis is determined by 



When the aerodynamic reference center coincides with the center of mass, the body 
aerodynamic moments are given by 



Riser and Suspension- Line Forces 

With the aerodynamic forces and moments known, the problem now is finding the 
forces and moments caused by the connecting lines. For this, consider a vector r^ 

(from the parachute to the payload as shown in fig. 1) which results in the following 
equations as developed in reference 8. 

= Rj + Pj - R 2 - P 2 ( 5 ) 
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is given by 



Figure 1. - Vector diagram of parawing 
configuration. 


The time rate of change of r^ 

dF l ^ ^ ^ ^ 4 . 

“dF - r l _ R 1 + P 1 " R 2 " p 2 ^ 

The vector p ^ is a vector fixed in the 

payload and can be expressed in the pay- 
load referenced frame as 

P 1 = X 1 i 1 +Y 1 3 1 + Z 1 A 1 (7) 

By taking the derivative of p^ with 
respect to time 

x ri + Y i'V z A “ "i x <>1 (8) 

Similarly 

P 2 ~ ^2 X P 2 (9) 


The time rate of change of and Rg is equal to the inertial velocity of points 1 
and 2. 



By combining equations (6), (8), and (10) 


r l =v l- v 2 + C °l Xp l' W 2 Xp 2 


( 10 ) 


(ID 
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The time rate of change of also can be expressed in the parachute reference frame 
as 


2* 1 2 + i 'ly, 2>2 +i 'lz,2 k 2 


( 12 ) 


By equating the expressions for r^ 


A — — 


‘'lx, 2*2 + ‘ly, 2 J 2 + ‘iz, 2 k 2 = V 1 " v 2 + "l X P 1 ’ "2 X p 2 


(13) 


The vectors and c x are already expressed in the axis system of the para- 
chute. Let [r] be a transformation matrix which transforms vectors from the axis 
system of body 1 to the axis system of body 2. (The equations for determining [r] 
are given in appendix A. ) Now 


‘lx, 2 S 2 + ‘ly, 2 J, 2 + ‘lz, 2^2 ' (" V x, 2 


“y,2 p 2z,2 + "z, 2 p 2y,2) 1 2 


(" v y, 

(- V z,2 
+ [ r il( V x, 


+ nPo™ 2)^2 
x, 2 H 2y, 2 T w y, 2^2x, 2^2 


2 z, 2^2x, 2 T x, 2^2z 

rtPc),, r> + t)P, 


+ Ci> 


1 ' “y, l^lz, 1 z, l^ly 
+ co_ ,p„ * - U) 


+ r i2( V y, 1 - ~z, rix, 1 


W z, l P ly, l) 
, l P lz, l) 


+ r 


13 ( V z, 1 + W x, l p ly, 1 " W y, l p lx, l) *2 


f 2l( V x, 1 + w y, l p lz, 1 ' w z, l P ly, l) 


+ r, 
+ r, 


22( V y, 1 


+ w z, l P lx, 1 W x, l p lz, l) 


23( V z, 1 + W x, l P ly, 1 " w y, l p lx, l)] h 


_ r 3!( V x, 


, + w .p, - w 1 p 1 

1 y, l^lz, 1 z, r ly, 1 ; 


+ r 


32 ( V y, 


+ 0 ) 


- Cl) 


+ r, 


33 


(L, 


1 1 w z, i M ix, 1 x, r iz 
1 + w x, l p ly, 1 " "y, l p lx 


,1) 

,i)] k : 


(14) 
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Equating the coefficients of 1 2 , j 2 , and & 2 results in the following three scalar equa- 
tions for the time rate of change of the components of each vector from the parachute 
to the payload. 


r lx, 2 V x, 2 W y, 2 P 2z, 2 + W z, 2 p 2y, 2 

+ r n( V x, 1 + W y, l P lz, 1 ' “z, l P ly, l) 

+ r i2( V y, 1 + w z, l P lx, 1 ” W x, l P lz, l) 

+ r q/V + a; , - co ,0. ,\ 

13V z» 1 x, l^ly, 1 y, l^lx, l) 

r ly, 2 = -V y, 2 _ W z, 2 P 2x, 2 + w x, 2 P 2z, 2 

+ r 2l( V x, 1 + W y, l P lz, 1 ' w z, l p ly, l) 

+ r 22( V y, 1 + w z, l p lx, 1 w x, l P lz, l) ^ 

23V z , 1 x, r ly, 1 y, rlx, 1) 

r lz, 2 = ” V z, 2 " w x, 2 P 2y, 2 + w y, 2 p 2x, 2 

+ r_,/V , + co ,p, .. - co ,p, A 

3 1 V X, 1 y, rlz, 1 z, r ly, l) 

+ r„ fv n + co ,p. 1 - co .p. A 

32V y> 1 z > 1 lx, 1 x, l'lz, 1) 

+ r„„/V , + co ,p, , - co .p, A 

33V z. 1 X, r ly, 1 y, l^lx, l) 


(15) 


The integration of these equations yields the displacement of each point on the payload 
with respect to the parachute reference frame. 

If the two points were connected by individual elastic lines such as r ^ in fig- 
ure 1, then the resulting forces transmitted to the parachute and payload are deter- 
mined by multiplying the spring constant of the line by its elongation. The dynamics 
of the system are then determined. The equations of motion of both the parachute and 
payload are integrated to yield the velocities and angular rates which are used to deter- 
mine the elongation of the connecting lines, thus giving the forces in these lines. 
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In most systems, however, the para- 
chute and payload are not connected by indi- 
vidual lines but by space frames composed 
of suspension lines with risers or harnesses 
connecting the confluence point of the sus- 
pension lines to the payload (fig. 2). In 
order to accurately obtain the system dy- 
namics, it is necessary to determine the 
force distribution in these space frames. 


Space Frame Forces 

To solve for these forces, consider 
the typical space frame as shown in fig- 
ure 3. Points 1 to 3 represent attach points 
on the payload, and the displacement of 
these points with respect to the parachute 
can be found by integrating the time rate of 
change of a vector from some arbitrary 
point on the parachute to the payload attach 
point. Points 4 to 8 are attach points of the 
suspension lines on the parachute. One of 
the basic assumptions of the analysis is that 
the canopy maintains a rigid shape such that 
the skirt of the canopy, where the suspen- 
sion lines attach, maintains a fairly con- 
stant shape. A review of flight test films 
indicates this to be a valid assumption for 
most canopy shapes over their normal 
angle -of- attack range. Small movements 
of the canopy attach point will have little 
effect on the analysis because of the long 
suspension- line lengths and the low spring 
constant of the suspension lines. Thus, the 
parachute attach points can be assumed to 
be fixed with respect to the parachute axis 
system. To reduce computations, the para- 
chute suspension lines are represented in 
the analysis by a reduced amount having an 
equivalent spring constant. The method 
developed, however, is applicable to any 
number of suspension lines and is limited 
only by the capacity of the computer used 
in the analysis. 

For a given frame, the forces in the 
members caused by a small displacement 
of the payload attach points can be deter- 
mined by one of the methods of indetermi- 
nant structures, if the frame members are 



Figure 2. - Suspension-line geometry 
of parawing configuration. 


A 



Figure 3. - Riser and suspension- line 
space frame. 
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considered to be linearly elastic. This approach assumes that, at any instant of time, 
the frame is in static equilibrium. This is a valid assumption as long as the frequency 
of oscillation between the parachute and payload is fairly low so that the acceleration 
of the frame does not prohibit its maintaining near-equilibrium conditions. The small- 
displacement theory of indeterminant structures is not valid, however, because the 
dynamics of the parachute payload system are such that the displacements of the 
payload attach points with respect to the parachute axis are not always small. 


Forces in Space Frames Subjected to 
Large Displacements 

An analytical technique is now developed which solves for the force distribution 
in the space frame subject to large displacements. This method follows that developed 
for small-displacement theory by Argyris (as discussed in ref. 9) but is expanded to 
incorporate large displacements. The following analysis considers a limited number 
of frame members, but can be readily expanded to incorporate any number of members. 


Let points 3 to 5 of figure 4 repre- 
sent points on the parachute which are 
fixed with respect to the parachute axis 

system ig, jg, kg. Points 1 and 2 rep- 
resent points on the payload which undergo 
displacements dx., dy^ and dz^ with re- 
spect to their initial position in the para- 
chute reference frame. The displacements 
of point o from its original position are 

dv d V and dz 0 - 

In the parachute reference frame, let 
Fx i’ Fy i’ 311(1 Fz i re P resent the three 

components of force at each movable point 
required to hold the frame in equilibrium 
when it is displaced from its unstressed 
position. Let a., /3., and y be the di- 
rection cosines of frame member i in the 
parachute reference frame and P^ be the 

axial force in member i. The forces Fx., 

i’ 

Fy^, and Fz. can then be written in terms 



Figure 4. - Space frame equilibrium 
forces. 
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of the axial forces and direction cosines as follows 




Figure 5. - Large displacement diagram. 


To solve for the force P. in terms of the 

large displacements, consider figure 5. A 
single member in planar motion is shown, 
but the theory readily extends to several 
members in three dimensional space. 

The axial force in member i is a 
function of its linear spring constant and 
elongation; thus 

p . = K.6. = Kjj[dx 2 (cos {. ) - dx^cos {,) 

+ dy 2 (cos - dy^cos ^/.J] 

- L. (1 - cos Ad) (17) 

_ l 0 
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In matrix form, equation (17) becomes 



The total angle change of frame member i from its unstressed condition is determined 
from 


12 2 
1 - cos A 6=2 ^ cos ^ - cos ^ + ^cos - cos ^ (19) 

By making this substitution and extending equation (17) to the three dimensional frame 
of figure 4 
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By substituting for P in equation (16), the force matrix can now be symbolically 
written as 


\F\ = [A] [K] 


[A] T id} - j^{L 0 }(p 2 +t 2 +a 2 )J 


(21) 


Thus, it is seen that the forces required to hold the space frame in equilibrium for a 
given displacement can be expressed in terms of the displacements of the points where 
the forces are applied and the original length and direction cosines of the frame mem- 
bers. 


Consider again the frame in figure 4. Points 1 and 2 are the payload attach 

points, and external forces are required at these points to displace the frame. 

Point o, however, is an interior point in the frame and no external forces are applied 

here. Thus, Fx , Fy , and Fz are virtual forces and vanish when the frame is in 
’ o o o 

static equilibrium. 

Performing the matrix multiplication required for finding Fx q , Fy Q , and Fz q 

yields the following three equations for a frame with m total members and n mem- 
bers from the payload to confluence point. 


■> 

Fx 0 = o - -£ (Vi 2 *^ + V A dy i + K iVi dz i) 

1=1 

+ £^[Vi 2dx o + VA dy o + W. + 2 Vi L i 0 ( p i 2 + T i 2 + a i 2 )] 

Fy o = 0 = + Vi 2dy i + Vi y i dz i) 

1=1 

m r- P -i / o 2 2\ 

+ £ + K A dy o + K A y i dz o + 2 K A L i n ( p i + T i + a i ) 

i=lL o 

Fz „ - 0 = -£( K i‘'i , 'i dx i + K i' 3 i y i dy i + K i y i 2dz i) 

1 = 1 ' ' 


> ( 22 ) 


rn 2 i / 2 2 2\ 

+ V K a y.dx o + K./3.r.dy o + + 5 K^L. Ip. + r, + J 

1 =1L U 
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Now assume that points 1 and 2 are given known displacements dx., dy., and dz^ The parameters a 

^i’ y V P i’ T i’ and °i are a11 functions of the variables dx Q , dy Q , and dz Q ; and equation (22) yields three 
equations of the form 



Thus, there are three equations in the three unknowns dx Q , dy Q , and dz Q . However, 

these equations are very nonlinear in the three unknowns, and their solution requires 
the use of an iterative method. The iterative substitution method of Gauss-Seidel was 
first attempted, but failed to converge for large displacements. 


Newton- Raphson Iterative Solution 

The Newton-Raphson method was next used and very good convergence for all dis- 
placements was obtained. Reference 10 explains the principles of the Newton-Raphson 
method, and its application to the present problem is explained below. Let the equa- 
tions for Fx q , Fy Q , and Fz q be denoted as the functions Fj, F 2 , and F g , or in 

matrix form as { F } . Define as <£ the 3 by 3 matrix composed of elements of the 
Jacobian 


a ( F r F 2 ’ F 3) 

a ( dx 0 i d y 0 , dz o) 


Then 


cj) = 


11 

*12 

*13 

21 

*22 

*23 

31 

*32 

*33 


(24) 


(25) 


The determinant of [<f>]. is the Jacobian of the system evaluated at step k. Equations 

for all F terms and all <fi terms are given in appendix B. The equations in appen- 
dix B were obtained by hand calculation and were later checked with a computer pro- 
gram using the FORMAC computer languages as described in reference 11. 

Now assume that initially at step k, dx Q , dy Q , and dz Q are known. If the pay- 
load now experiences an additional displacement, the new values of dx Q , dy Q , an d dz Q 
can be computed as follows. 

Let 


/ 5dx 


dx \ 

°l 


0 l 

5dy o J 

= corrections to previous values of 


[ 6dz o 


dz o) 


(26) 
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Then by the Newton- Raphson method 


where 


dx 

° 


dx 

o 

| 

6dx 

o 

dy 

J o 

■ 

dy o 

+ 

5dy o 

dz 

o 

Ik- , 1 

dz 

o 

. 1 

6dz 

o 


= W k+ 1 


6dx 

o 

6 d y 0 


= -Wk"'f p tk 


6dz 

o 


k 


(27) 


(28) 


The matrices [0]^ 1 and |F} k are evaluated using the values of 


Using the new values of |d(, the process can be repeated until the solution con- 
verges as evidenced by the residual of the functions approaching zero. 

The Newton- Raphson method converges rapidly, usually requiring only three or 
four iterations for a typical 10- to 20-member frame. The initial approximation for 
|d} k is taken as the unstressed position of the junction, and subsequent approximations 

are taken from the position at the end of the previous integration step. In general, 
these are good approximations and are a factor in the rapid convergence of the Newton- 
Raphson method. The primary problem is determining the partial derivatives of { Ff, 
but once they are found, the equations can be adapted to handle any number of frame 
members as long as only one junction exists per frame. Equations for the partial de- 
rivatives are given in appendix B and are valid for any frame so long as only one con- 
fluence point exists in the frame. 

After the displacements of point o are found, the forces in each frame member 
can be determined by the matrix operation of equation (20). 


dx \ 

o I 

dy Q j at step k. 
dz o) 



Since the direction cosines of each frame member in the parachute axis system 
have been computed, it is possible to determine the force components of each member 
in this system. The force components of the frame members attached to the payload 

T 

can be found by multiplying by the transformation matrix [T] . If these force compo- 

nents are now multiplied by the moment arms from the parachute and payload centers 
of mass, the moment contributions caused by the frames can be calculated. These 
forces and moments can now be added to the aerodynamic forces and moments to cal- 
culate the six-degree-of -freedom motion of both the payload and parachute. 


CONCLUDING REMARKS 


The mathematical model as presented will permit a more detailed analysis of the 
flight dynamics of general parachute-payload systems than has been possible in the 
past. A particular advantage of the two-body approach over a single rigid-body analy- 
sis is that the aerodynamic damping coefficients are much less critical because of the 
large distance between the canopy center of gravity and the combined canopy -payload 
center of gravities. In the single-body analysis, this distance is combined with the 
angular velocity of the system to give the induced velocity term necessary for damping 
equations. In the present analysis, damping coefficients have less effect because the 
induced velocity of the canopy is calculated in the equation of motion. Since damping 
coefficients are by far the most difficult to determine in wind tunnel testing, the advan- 
tage of using this method is obvious. In addition, the two -body approach allows relative 
motion between the parachute and payload which can result in significant changes in the 
inertia matrix of the equivalent single rigid body. The single-body approach does not 
account for these inertia changes. 

The mathematical model as presented has been programed in FORTRAN IV and V 
languages for the IBM 7094 and Univac 1108 computers. The qualitative results ap- 
peared very good, but lack of experimental data prevents a complete evaluation of the 
program. Although studies are currently underway to determine aerodynamic and mass 
properties of lifting parachutes (refs. 12 to 15), as yet, they have not been defined well 
enough to have complete confidence for use in six-degree-of-freedom simulations. A 
program listing of the model will be supplied on request. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
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APPENDIX A 


EQUATIONS FOR [P] 


The transformation matrix [T] which transforms vectors from the axis system 
of body 1 to the axis system of body 2 is determined from 


(cos 9 cos \p) 


(cos 9 sin ip) 


(-sin 9) 


[r] = 


(-sin xp cos 4> _ (cos cos x p _ (sin cos 9) 

+ sin <p sin F cos xp) + sin $> sin 9 sin xp) 

(sin xp sin (p_ _ (-sin (p cos xp _ (cos cp cos F) 

+ cos sin F cos i p) + cos ?> sin 9 sin xp) 


The angles <p, 9, and xp are the relative^ Euler angles of body 2 with respect to 
body 1 applied in the order of rotation xp,9,(p. These relative Euler angles are deter- 
mined by integrating the equations 

^ - K 2 - “x, l) + *“ * [("y, 2 - "y, l) Sln ? + K, 2 ’ “z, l) cos ^ <A2> 

^ = ("y, 2 - "y, l) C ° S ? ' (" 2 , 2 ' “z, l) sto ? (A3) 

$ = [("y, 2 ■ "y, l) s “ ? + ("z, 2 • “z, l) cos ^ * cos ? < A4 > 


CO , 
X, 1 


u> 1 
X, 1 

co , 

y, 1 

= [r] 

CO , 

y> 1 

<x> 1 

z, 1 


co 1 
z, 1 


APPENDIX B 


EQUATIONS FOR EVALUATING ELEMENTS OF THE JACOBIAN 


" /k. AX.V ' ^AX.AY.dy. K.AX^Z.dz A 

1 ~2-j\ T 2 + T, 2 L. 2 I 


. 1 \ L. 

1=1 \ l 


L. 
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E 

i=l 


K. AX.L. 


K.AX, dx Q K. AX.AY.dy 0 K.AX.AZ.dz,, 1 1 ‘ 0 /. 2 . . 2 . , 2\ 

— T ~l + T2 + T~2 Ifi +T i +a i ) 


L. 

1 


L. 

l 


L. 

l 


(Bl) 


F = - 
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iL/K.AX.AY.dx. K.AY. 2 dy. K.AY.AZ dz\ 
V M 1111 111,1111' 


i=l \ L i 
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L. 
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i=l 


K. AY.L. 
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L. 

l 
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L. 
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Jl* /2K.AX.dx. 2K.AX. 3 dx. K.AY.dy. 2K.AX. 2 dy.AY. 
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